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ABSTRACT 

We present a Spitzer Space Telescope imaging survey of the most massive Galactic 
globular cluster, uj Centauri, and investigate stellar mass loss at low metallicity and the 
intracluster medium (ICM). The survey covers approximately 3.2 x the cluster half-mass 
radius at 3.6, 4.5, 5.8, 8, and 24 fim, resulting in a catalog of over 40,000 point-sources 
in the cluster. Approximately 140 cluster members ranging 1.5 dex in metallicity show 
a red excess at 24 /im, indicative of circumstellar dust. If all of the dusty sources are 
experiencing mass loss, the cumulative rate of loss is estimated at 2.9 - 4.2 x 10"'^ Mq 
yr~^, 63% - 66% of which is supplied by three asymptotic giant branch stars at the tip 
of the Red Giant Branch (RGB). There is little evidence for strong mass loss lower on 
the RGB. If this material had remained in the cluster center, its dust component (> 1 
X 10~^ -^0) would be detectable in our 24 and 70 fim images. While no dust cloud 
located at the center of wCen is apparent, we do see four regions of very faint, diffuse 
emission beyond two half-mass radii at 24 iim. It is unclear whether these dust clouds 
are foreground emission or are associated with u) Cen. In the latter case, these clouds 
may be the ICM in the process of escaping from the cluster. 

Subject headings: globular clusters: individual (a; Centauri) — infrared: general — stars: 
AGB and post-AGB — stars: mass loss — stars: variables — ISM: evolution 



Introduction 



Despite their old age and low metal content. Galactic globular clusters (GCs) are now known 



to be dust factories (Ramdani Sz Jorissen 2001 Origlia et al. 2002 Evans et al.|2003 Boyer et al 



2006 van Loon et al.||2006a, Lebzelter et al.|2006 Ita et al^|2007 |. This has important implications 
for our understanding of mass loss from low-mass and metal-poor stars and for the replenishment of 
the interstellar medium (ISM) within old, metal-poor stellar systems, such as GCs, dwarf spheroidal 
(dSph) galaxies, and galactic halos. 

All stars more massive than ~ 0.8 Mq lose a significant fraction of their mass before they 
deflagrate or leave a compact remnant (see reviews by |Iben Renzini|1983]|Chiosi k Maeder|1986 |). 
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In the final stages of evolution, most stars enter higher stages of nuclear burning, and the products 
of this nucleosynthesis are transported to the stellar surface through convection or rotation. The 
chemically enriched mass that is recycled (with some delay) into the ISM has a significant impact 
on any further star formation, and drives the chemical evolution within galaxies. In particular, 
the enrichment of the ISM by dust is of great importance, as grains play important roles in many 
ISM processes, including the formation of molecular hydrogen ( Gould Sz Salpeter||1963 1 and planet 
formation. The winds of low and intermediate-mass stars on the upper red giant and asymptotic 



giant branches (RGB and AGB) and massive red supergiants are prolific dust producers (Gehrz 



19891. The dust is an integral part of the driving mechanism of the most massive winds (e.g., Gehrz 



Sz Woolf 1971 Jura Sz Morris 1985 Gail Sz Sedlmayr 19871, and our understanding of red giant 



mass loss relies heavily on the details of the dust condensation process (e.g., Salpeter||1974 Ganger 



et al. 


1990 


Hofner &; Dorfi 


1997 


(for a review, see 


van Loon 


2006). 



which is poorly understood in particular at low metallicity 



GCs are exquisite laboratories in which to study stellar winds because they represent the 
closest match to a single stellar population at metallicities that span over two orders of magnitude 
down to < 1% solar (cf. Gratton et al.|[2"004 |. Mass loss remains important throughout the history 
of a GC (van Loon &z McDonald 20081; even in old GCs, stars lose ~30% of their mass, most or 
all of which is removed from the cluster via a host of mechanisms including ram-pressure stripping 
during Galactic plane crossings. RGB mass loss may be responsible for the blue horizontal branches 



observed in metal-poor GCs (Rood 19731, and evidence for mass loss is seen in the line profiles of 



chromospherically active or pulsating red giants in GCs (Dupree et al. 1984; Cacciari et al. 2004 



McDonald k van Loon 2007 Meszaros et al. 2008). 



There is mounting evidence for mass loss from GC red giants in the form of dusty winds 



(Ramdani Sz Jorissen 2001 Origlia et al. 2002). The first spectroscopic identifications of silicate 



dust are in the tip- AGB star VI in the massive GC 47Tucanae (van Loon et al. 2006a I and 
subsequent identification of minerals in other AGB stars in 47 Tuc with Spitzer IRS observations 
(Lebzelter et al. 2006). With [Fe/H] = —0.7, 47Tuc is a metal-rich GC, providing an environment 
similar to that found, for instance, in the Small Magellanic Cloud. To probe uncharted terrain, it 
is important to study dust in GCs at lower metallicity. The surprising detection of both interstellar 
and circumstellar dust and gas in the metal-poor GC M15 dEvans et al.]2003| [B^ et al.|2006 



van 



Loon et al.||2006b | demonstrates that dust forms even in extremely metal-deficient environments. 



1.1. uj Centauri 



The most massive Galactic GC, u; Centauri, contains enough stars in the dusty wind phase 
for an empirical reconstruction of the evolution of mass loss. Its gravitational well (wcsc = 44 
km s~^; Gnedin et al. 2002) could retain interstellar matter if winds are slow. wCen is located 



approximately 15° out of the Galactic plane at a distance of ~5 kpc (Lub 2002 van de Ven et al 



2006[ [Del Principe et ar]|2006[ ), so it is possible to probe the stellar population to the cluster core. 

The bulk of the stars in uj Cen have [Fe/H] = -1.7 ( |Norris|1996| [Smith et al.|2000[ ), a metallicity 
much lower than the environments in which circumstellar dust has traditionally been studied. 
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Complications arise because u) Cen also harbors a small fraction of stars up to an order of magnitude 
more metal-rich ( Lee et al.|1999 Pancino et al.|20"00 20021 and an intermediate population of stars 
with an abnormally high helium content ( Norris|1996 2004 1. Red giants in uj Cen display a range in 
surface abundances, from oxygen-rich, M-type stars with titanium oxides, stars enhanced in CH or 
in cyanide (CN), to genuine carbon stars with molecular carbon (C2; van Loon et al.|[2007| ). These 
disparate populations cause difficulties in stellar evolution studies of uj Cen, but since stars from the 
metal-rich and metal-poor populations can be distinguished by means of optical spectroscopy, lo Cen 
presents us with a unique opportunity to study the dust and mass loss within a single environment 
as a function of atmospheric chemistry. 

Consensus has yet to be reached regarding the formation history of the different subpopulations 



(Sollima et al. 2005 Stanford et al. 2006; Villanova et al. 2007). The spread and peculiarities in 



elemental abundances in uj Cen have been interpreted as evidence that this GC may be the remnant 



nucleus of a tidally disrupted dSph (Zinnecker et al. 1988 Freeman 19931. If so, w Cen may offer 



us insight into the mass loss and chemical enrichment in the much more distant nucleated dSphs, 
and perhaps dSphs in general. 



Here, we present NASA Spitzer Space Telescope ( Werner et al.||2004 Gehrz et al.||2007 | images 
of u> Cen and a catalog including the magnitudes of over 40,000 point-sources from 3.6 to 70 fim. The 
24 fim mosaic, which enables reliable detection of individual circumstellar dust shells near the peak 
wavelength of their spectral energy distributions (SEDs), provides the most dramatic improvement 



over previous data. We make use of proper motion (van Leeuwen et al. 2000 1 and radial velocity (van 



Loon et al.|2007 | measurements to separate cluster members from the substantial foreground stellar 
and background galaxy populations. Optical spectroscopic information permits us to compare 
the contributions from stars of different composition and evolutionary status to their infrared 
(IR) emission and dust production. Our observations are discussed in §[2] photometry and the 
a; Cen point-source catalog are described in §[3j The color- magnitude diagrams (CMDs), stellar 
populations, mass loss, and the intracluster medium (ICM) dust are discussed in §|4j and our 
conclusions are presented in §[5j 



Observations 



We used the NASA Spitzer Space Telescope to obtain image maps of the GC uj Cen (Program 
Identification [PID] 20648, PI: van Loon) with the Infrared Array Camera (IRAC; [Fazio et al.|2004 | 



and the Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004 1 . Observations are 
centered at R.A. = 13'"26"^45';9, decl. = -47°28'37'.'0 [J2000.0] and map the cluster to at least 3.2 
times the half-mass radius (rhm = 4.8'; Harris 19961 at each wavelength. Observation details are 
summarized in Table [T] 

Images centered at 3.6, 4.5, 5.8, and 8 /im were obtained with IRAC on 2006 March 26 UT. 
We observed in high dynamic range mode, enabling us to correct saturated sources in the long 
exposures (10.4 s) with their counterparts in the short exposures (0.4 s) while still achieving a high 
signal-to-noise ratio for the fainter sources. Images were taken at nine dither positions in a cycling 
pattern to build redundancy against outliers and artefacts. The resulting image map contains a 
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total of 441 images and covers an area approximately 35' x 35' around the cluster center. 

MIPS images at 24 and 70 fim were obtained on 2006 February 22 UT and 2006 March 02 
UT with seven scan legs of approximately 0.5° in length. The resulting mosaics combine 4,620 
images for each wavelength and cover a total area of approximately 35' x 55'. We designed our 
MIPS observing campaign with the intention of obtaining only 24 /im data, and the MIPS 70 and 
160 /um data were obtained automatically (see Rieke et al.||2004 |. Due to the non-functioning half 



of the 70 nm array, the 70 /umdata do not have complete coverage. The warm state of the telescope 
during the observations saturated the 160 /im array, rendering the 160 /xm data useless. MIPS 
observation details are summarized in Table [TJ 

Raw MIPS and IRAC data were processed with the Spitzer Science Center (SSC) pipeline, 
version S13.2.C|5 The Basic Calibrated Data (BCD) were post-processed with portions of the 2006 



March 01 version of the SSC Legacy MOPEX software ( [Makovoz k Marleau|[2005[ ). For IRAC, 
after applying an array distortion correction, we implemented the MOPEX background-matching 
routine to minimize pixel offsets in overlapping areas of the mosaics and the MOPEX mosaicker 
to eliminate cosmic rays and other outliers. No further cosmetic corrections were applied to the 
BCDs before performing point-spread-function (PSF) photometry. The final mosaics have a pixel 
size of 1'.'22 pix"^. A three-color image of the first three IRAC bands is presented in Figure [l^, in 
which very few sources appear with a red or blue excess. 

For the MIPS 24 /um data, we created and applied scan-mirror-dependent and scan-mirror- 
independent flat fields for all images excluding the bias boost images, any image within 10' of the 
cluster center, and any image containing a very bright source. A 2-degree polynomial fit was applied 
over the entire image to subtract the zodiacal background. At 70 fim, we first performed pixel time- 
filtering and column filtering with SSC contributed software from D. Fadda and D. Frayei[^ We 
then implemented the MOPEX background-matching and mosaic routines at both 24 and 70 /im 
before proceeding with point-source extraction. An IRAC and MIPS three-color image is presented 
in Figure [l]3, in which several cool dusty objects appear green at 8 fim and red at 24 fim. 



3. Photometry 

3.1. Point- Source Extraction 

Point-source extraction for 3.6 - 24 /im was done with the DAOphot II photometry package 
. PSFs were created using a minimum of 20 isolated stars. Sources brighter than 4 a 
over the background were chosen for extraction from the mosaic image. Any remaining outliers and 
other extended objects were eliminated from the sample with a sharpness cut-off, which identifies 
sources broader or narrower than the PSF. Photometric completeness at each wavelength was 
determined through false star tests; 50% and 90% completeness limits are hsted in Table |2j 



(Stetson! 1987 



^http:/ /ssc. spitzer. caltech.edu/archanaly/plhistory/ 
^http:/ /ssc. spitzer. caltech.edu/mips/filter/ 
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The final fluxes and Hux uncertainties are color-corrected according to the IRAC and MIPS 
Data Handbooks versions 3.0 and 3.2.1, respectively {Spitzer Science Center 2006a|b ), using a 
5000 K blackbody, which is the approximate temperature of a typical RGB star. Fluxes in the 
point-source catalog described in § 3.2 have not been corrected for reddening. Where the reddening 



correction has been applied for analysis of the data, Ax values are taken from Indebetouw et 
al. (20051 and E{B — V) = 0.11 (Lub 2002|. The pixel-phase-dependent photometric correction 



described by Reach et al. (20051 was also applied to 3.6 fj,m fluxes and their uncertainties. This 
correction was applied only to 3.6 fj,m because it is not well defined for the other ^piteer wavelengths. 
The flux uncertainties that we report in the Spitzer catalog described in section §|3.2|include both the 



uncertainties quoted by DAOphot and the calibration uncertainties specified by Reach et al. (20051 
for IRAC and Engelbracht et al. (20071 for 24 iim. All fluxes were converted to stellar magnitudes 
using the zero fluxes reported in the Spitzer Data Handbooks. The photometric uncertainty as a 
function of source magnitude is given in Figure |2} 



3.2. The Point-Source Catalog 



We have created a point-source catalog of > 40,000 wCen objects detected at 3.6, 4.5, 5.8, 8, 
and 24 /um. An extraction from this catalog is given in Table |3j where column one is the source 
number, column two is the source ID, columns 3-7 are Spitzer magnitudes, column 8 is the 



membership probability (van Leeuwen et al. 20001, columns 9 and 10 list source type labels, and 



column 11 lists the Leiden identifiers (LEID) from van Leeuwen et al. (20001. The source ID follows 



the standard Spitzer naming convention, giving the truncated (J2000) coordinates. The complete 
Spitzer catalog is available electronically through the Centre Donnees de Strasbourg (CDS). 

Cluster members and non-members were identified by cross-correlating the Spitzer source co- 



ordinates with those in the van Leeuwen et al. (2000) proper motion catalog. The van Leeuwen 



et al. (20001 catalog is not as deep as our survey nor is the coverage as wide, so many Spitzer 



sources do not have proper motion measurements, van Leeuwen et al. (20001 derived positions and 



proper motions by cross-referencing with the International Celestial Reference Frame as defined by 
the Hipparcos catalog and extended to a denser grid of stars by the Tycho-2 catalog. The coordi- 



nates of the Spitzer sources are systematically offset from the van Leeuwen et al. (20001 sources, 
with mean offsets of Aa = -1-0.40" and A6 = -1-0.25" and standard deviations of 0.15" and 0.10", 
respectively. Coordinates in the Spitzer catalog are not corrected for this offset. 



If a source in any IRAC band was not associated with another IRAC, MIPS, or van Leeuwen 



et al. (2000) object, we considered it a false detection and dropped it from the Spitzer catalog. 



MIPS 24 fim sources that do not have IRAC or optical counterparts remain in the catalog since the 
very red SEDs of many MIPS sources peak in the mid-IR and are thus fainter than the sensitivity 
limits in the bluer Spitzer bandpasses. 

Also included in the Spitzer catalog are notes in columns 9 and 10 indicating if the source is 
a galaxy, AGB, anomalous RGB (RGB-a), or Horizontal Branch (HB) star, based on its locations 



on the CMDs (see §4.1). The catalog also indicates sources that are super-Ba- or CN-rich and/or 



are M-type (M), carbon (C), or post- AGB stars (P), as identified by van Loon et al. (2007). 
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The differing resolution between IRAC and MIPS has caused some sources that are isolated 
in IRAC to be blended in the 24 mosaic, which produces an artificially red [8] - [24] color. To 
determine which sources are potentially affected by blending, we convolved the 8 //m mosaic with 
the 24 ^m PSF and identified sources with elevated fluxes in the convolved mosaic compared to 
the unaltered mosaic (Fig. [s]). These sources are marked with "B24", for "blended at 24 nia" , in 
column 10 of the Spitzer catalog. 



3.3. Luminosity Functions and Completeness 

Luminosity functions for epoch two are displayed in Figure |4j The distance modulus used to 
determine absolute magnitudes is 13.58, given a distance of 5.2 kpc ( [Del Principe et al. 20061. In 



the top panels, known non-members are excluded and of the remaining sources, only those that are 
in at least two adjacent wavebands are included. The IRAC luminosity functions show a strong 
change in slope near an absolute magnitude of -2 that may be the RGB bump or the RGB clump, 
both of which are located at absolute magnitudes of ss— 2 - in the near-IR data of several GCs 
( Valenti et al.|[2007 ), depending on the metallicity of the cluster. The HB is also visible at 3.6, 4.5, 



and 5.8 fxm just beyond an absolute magnitude of zero as either a bump or change in slope. 

The tip of the RGB (tRGB) is difficult to determine from the luminosity functions, as the 
source counts are low at absolute magnitudes brighter than —5. The sparse population of AGB 
stars above the tRGB is expected in GCs since this phase of stellar evolution is particularly short- 
lived and also because the very low mass of stars in GCs (Magb — 0.6 — 0.8 Mq) prevents 
AGB luminosity from climbing much above the tRGB. Nevertheless, the source counts of known 
cluster members (lower panels of Fig. |4]) do drop by a factor of 3 to 4 past the bins, which are 
0.2 magnitudes wide, centered at absolute magnitudes of —5.9, —6.1, —6.1, —6.1, and —6.3 at 3.6, 
4.5, 5.8, 8, and 24 //m, respectively. The K-band tRGB for a; Cen is located at ~— 6.04 it 0.16 



magnitudes (Bellazzini et al. 2004). At 3.6 fiui, the tRGB location is —6 magnitudes for many 



stellar populations, including the Large Magellanic Cloud (LMC; van Loon et al.|[2005 ) and smaller 



stellar populations such as the Local Group dwarf irregular galaxy WLM ( Jackson et al.|[2007 |. 



Photometric completeness in each waveband was determined by performing false star tests, 
for which 2% of the total stellar population that mimicked the luminosity function was added to 
each mosaic at random positions and PSF photometry was performed to attempt to recover these 
false stars. The test was carried out 100 times for each wavelength. Dotted lines mark the 50% 
completeness limits in the upper panels of Figure |4j which are listed in Table [2} The dashed line in 
Figure |4] marks the sensitivity limit at each bandpass, and shows that the photometry is the most 
incomplete at 3.6 and 4.5 fim, which is likely due to severe stellar crowding in the cluster core at 
these wavelengths. Photometry at 5.8 - 24 /im is 50% complete to within two magnitudes of the 3 
a sensitivity limits. 
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4. Discussion 



4.1. The Color-Magnitude Diagrams 



Figure [s] displays the optical ( van Leeuwen et al.||2000" |, the near- IR dSoUima et aT][2004l ), and 
our mid-IR CMDs. The CMDs are corrected for reddening, as discussed in §|3] with Ab values from 



Rieke &z. Lebofsky (19851, but the effect of reddening on the mid-IR magnitudes is minimal. The 



IRAC and MIPS CMDs shows very little structure aside from the RGB, and at longer wavelengths, 
a very broad, faint red plume. Mid-IR isochrones that include all known RGB/AGB properties, 
molecular bands, and dust considerations have yet to be developed. For reference, the brightest 
IRAC sources are labeled with large filled circles in Figure [5j and the three brightest M-type stars 
are also labeled with their LEID numbers. Bright IRAC sources that are not labeled are not cluster 
members. Figure [6] shows the same CMDs with different source types labeled in color. Black points 
are confirmed cluster members (confirmed non-members are not plotted) and gray points, which 
are mostly located at the faint end of the RGB, do not have matches in the proper motion catalog. 

Absorption in molecular bands may be responsible for some of the spread towards red colors 
on the RGB, but cannot cause colors as red as those of the dustiest stars. However, while many 
of the very red stars on the mid- to lower- RGB in Figure |5] may be dusty, it is very likely that 
some of the these sources are red due either to mismatches across bands, photometric errors, or 
contamination from cluster non-members. 

The B - [3.6] CMD clearly displays the RGB, AGB, and HB. Also visible as a sparsely pop- 
ulated branch to the red of the RGB is the RGB-a (magenta points). The plume of stars that 
falls intermediate to the RGB and the HB is mostly due to foreground contamination ( |Bellazzi"iii 
et al. 20041, however the brighest cluster members in this plume (-B < 13, 1 > B - [3.6] > 3) 
include post-AGB stars such as LEID 16018, 30020, 32015, 32029, and 43105. Note that while 



van Leeuwen et al. (20001 found that 16018 and 32029 are cluster non- members, van Loon et al 



(20071 determined through optical spectroscopy that they do in fact belong to a; Cen. As such. 



their membership probabilities have been set to 100% in the Spitzer catalog. 



We identified HB stars (blue points) by their positions on the optical CMD (van Leeuwen et 
al] |2000[ ). At IR wavelengths, the HB becomes indistinct from the RGB. The stellar density profile 
shown in Figure [7] shows a paucity of HB stars in the center of the cluster. The non-central peak 
of the HB density profile is explained by the dominance of the cluster luminosity by bright AGB 
and RGB stars in the most crowded regions. 

RGB-a stars are those located on the red branch of the RGB in optical and near-IR CMDs 



(see Fig. 5 of Sollima et al. 2004). The offset of the RGB-a from the main RGB in optical and 
near-IR CMDs has been attributed to a metallicity difference between the RGB-a and the regular 
RGB ( |Lee et al.||1999t|Pancino et al.||2000||2002t|OrigUa et al.||2003[ ). Except for a slight blue excess 
in the [3.6] - [4.5] colors, no metallicity separation is evident in the Spitzer CMDs. 



- 8 - 



4.2. The Asymptotic and Red Giant Branches 



Inspection of the luminosity functions (Fig.|4]) suggests that the stehar count begins to decrease 
at an accelerated rate near an absolute magnitude of —5 at each wavelength. We chose this 
magnitude as the divide at all Spitzer wavelengths between mass-losing AGB candidates and the 
rest of the stellar population. We identify 75 mass-losing AGB candidates, several of which have a 



red excess, although mass-loss can occur even without a dusty circumstellar envelope (McDonald Sz 



van Looii||2007 ). The AGB population is slightly more centrally concentrated than the total stellar 
population (Fig. [7| , which is probably due to decreasing completeness of fainter sources towards 
the core. 

The locations on the CMD of stars that are identified as CN- and super-Ba-rich by optical 
spectroscopy ( van Loon et al.|2007 | are marked in Figure [Sji. The term super-Ba-rich distinguishes 
stars with enhanced Ba through pre-enrichment from those that are enriched beyond this, probably 
due to internal mixing processes. The super-Ba-rich stars dominate the stars brighter than [24] ~ 



8.5, supporting the suggestion by van Loon et al. (2007) that these stars are enriched in Barium 



due to third-dredge up, and therefore delineate the thermal pulsating AGB. The super-Ba-rich 
stars may also be enriched in carbon, but not so much as to warrant their classification as genuine 
carbon stars. 

Figure [9] shows the [8] - [24] CMD for cluster members, non-members, and sources without 
proper motion measurements. Sources affected by blending at 24 /^m are not plotted. The location 
on the CMD of the majority of sources without proper motion measurements suggests that most 
are likely galaxies. Stars brighter than [24] ~ 10 with [8] - [24] > 0.1 likely harbor dust, while 
sources fainter than [24] ~ 10 are dusty if [8] - [24] > 0.2. Also marked in Figure |9] and listed 
m Table |4] are 16 stars of interest, identified by [van Loon et aT] ( |2007[ ) as six M-type stars (open 
squares), five carbon stars that show C2 bands (open triangles), and five post- AGB stars (closed 
circles), including Fehrenbach's star (Fehrenbach &; Duflot 1962) near [24] = 9.25. 



van Loon et al. (2007) obtained optical spectra for five of the extremely red faint cluster 
members ([24] > 10 and [8] ~ [24] > 0.8); these stars range from [Fe/H] = -1.5 to [Fe/H] = 
—2, and two of the stars are abnormally hot (T ~ 7000 K). The remaining three stars appear to 
be normal RGB stars. Visual inspection suggests that none of these five sources are affected by 
blending at 24 fim, however, it is possible that HB stars or RGB stars that have little or no emission 
at 24 fim are spatially coincident with background galaxies, causing them to appear to have red 
colors. This may also be the case for the other «10 faint, extremely red stars that are matched to 
cluster members. 

Spitzer studies of GCs that do not have proper motion or radial velocity unambiguous con- 
firmation of membership will suffer from contamination from red field stars mistaken for dusty 
cluster stars. Many of the bright 24 //m sources in our images are field stars, including a variable 
M-type star with a period of 509 days (LEID 34041). Proper motion measurements allowed us to 
identify almost all bright and red cluster members in Figure [9] against field stars, with only one 
bright star remaining for which no proper motion information is available. This bright source is 
most likely LEID 55017, which is a confirmed cluster non-member. The proper motion of LEID 
55017 is large enough that the Spitzer coordinates are ~4" different from the van Leeuwen et al 
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(20001 coordinates. 

It is expected that metal-rich stars form more dust in their winds than their low metallicity 
counterparts due to the higher abundance of condensable material. We have metallicity information 
from van Loon et al. (2007) for 32 of the «140 dusty stars in Figure |9j The brightest potentially 
dusty stars in the cluster range —2.25 < [Fe/H] < —1.25, i.e. nearly the entire breadth of metallicities 
encountered in uj Cen. In fact, of this small sample, the reddest stars have [Fe/H] < —1.5 (Fig. 10) 
and two of the three brightest dusty stars (LEID 33062, 44262, 35250) have [Fe/H] < -2.0, which 
suggests that dust production is not inhibited at these low metallicities. 



4.3. Mass-Loss Rates 

Since mass-loss rates (M) scale directly with dust optical depth (r), all dusty stars have the 
potential to drive a wind ( Ivezic &: Elitzur|l995 ). However, the value for M derived from the optical 
depth also scales with the square root of the luminosity, so the brightest stars contribute a larger 
fraction of the cluster's cumulative mass-loss rate than the fainter stars. Of the ~140 red cluster 
members that remain after eliminating blends and non-members such as LEID 34041, which, as 
a bright M-type star with a period of 509 days, would dominate the cumulative mass-loss rate of 
the cluster if included. We used the AGB models computed by Groenewegen ( 2006^ to determine 
the mass-loss rates corresponding to the color of the red member stars with 0.15 > [8] - [24] > 



3.1. Luminosities were determined by scaling to the Groenewegen (2006) mass-loss tracks for AGB 



stars with temperatures of 3297 and 3850 K and with silicate and AlOx dust compositions, however 
differing bolometric corrections for stars with different temperatures and metallicities are not taken 
into account. For reference, the colors corresponding to two sample optical depths are marked on 
the CMD in Figure gand listed in Table [HJ The final mass-loss rates (Mpinai) were scaled according 
to the following relationship: 



Ml 



Mg J 



Final 



10 km s 



-1 



0.005 



L 



UoOO Lp 



■0 



(1) 



where Mq is the unsealed mass-loss rate from Groenewegen (2006), v is the wind velocity, -i/^ is the 
dust-to-gas mass ratio ~ 0.005 at solar metallicity), and L is the luminosity. Following 



van 



Loon (2000) and Marshall et al. (2004), the velocity scales approximately as: 



10 km s 



-1 



0.005 



L 



10,000 Lq 



(2) 



From equation (2), the wind velocities of even the strongest mass-losing stars in a; Cen are slow 
(f < 2 km s^"*^). Winds this slow may have difficulty driving a massive wind; it is possible that 
the wind is driven by some other mechanism and that the dust is a by-product. McDonald &: van 
Loon (2007) suggest wind velocities of 5 - 10 km s^^ for AGB stars in GCs; if wind velocities are 



5x higher, the mass-loss rates would increase by the same factor. 
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The Groenewegen (2006) models were calculated assuming spherical circumstellar dust shells, 



a constant mass-loss rate, and photospheres with solar metallicity. The last of these assumptions 
affects the stellar colors, which we used to determine the mass-loss rates. While we do take metal- 
licity into account in the dust-to-gas ratio, increasing the mass-loss rates by a factor of seven, the 



metallicity dependence is not yet well established. The mass-loss rates computed using the Groe- 



newegen (20061 models should therefore be considered a rough estimate. Bolometric luminosities 



of the uj Cen sources can be estimated through fitting the SED or spectrum of each star, facilitating 
a better estimate of the real mass-loss rates (I. McDonald et al., in preparation). 



Assuming that all stars have the metallicity that comprises the majority of the population 

~ 10~^), we find that 



and that ip scales with metallicity (van Loon 



2000 



2006 



[Fe/H] 



-1.7, ^P 



the three brightest M-type stars (Fig. [9|) contribute a total of M(gas -|- dust) = 1.8 - 2.8 x 10""^ 
Mq yr~^, depending on the temperature and dust composition. The combined mass-loss rate of 
these three stars is an order of magnitude less than the mass-loss rate measured for VI in 47 Tuc 



with IR spectroscopy (van Loon et al. 2006a I, suggesting either that we have underestimated the 

and 
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mass-loss rates or that the variance of M with metallicity is not a linear relation. Figure 
Table |6] shows that the cumulative mass-loss rate (S M) of the ~140 dusty stars together is less 
than twice the combined mass-loss rate of the three brightest M-type stars (LEID 33062, 44262, 
35250). However, the fainter stars may be contributing more to the total mass-loss rate of the 
cluster if their dust-to-gas ratios are smaller. 

To compare S M to the number of dusty stars present at each magnitude, we also plot the 



percentage of stars that are dusty (/d) in Figure 11 , binned by half magnitudes. Jd remains roughly 
constant at 10% - 20% for [24] > 9.5, causing the cumulative mass-loss rate to slowly increase at 
[24] > 8. The increase in S M at [24] ~ 11 is due to the very red sources ([8] - [24] > 1.0) at these 
magnitudes. It is possible that many of the faint stars with red excess are not truly dusty (see 



^4.2). Since the cumulative mass-loss rate does not increase much over more than 2 magnitudes 



on the upper-RGB, it is unlikely that even fainter stars, which are generally warmer and do not 
pulsate strongly, contribute significantly to the dust-traced mass loss. 



In a Spitzer study of 47 Tuc, Origlia et al. (20071 find that red stars occupy the CMD from the 



tRGB down to the HB and conclude that mass-loss is occurring at a significant rate along the entire 
RGB. However, at least 25% of the red sources in a; Cen have proven to have red colors because 



of blending effects at 24 ^um (see § 3.2). Figure 11 shows that these blended sources can increase 
the cumulative mass-loss rate by up to 17%. In more crowded clusters like 47 Tuc, blending will 



affect an even larger percentage of red stars. Indeed, it is evident in Figure 2 of Origlia et al. (2007) 



that the CMD covering the central region of the cluster is affected by blending, since the many red 
sources claimed to exist between -2 < Mboiometric < have no equivalent in the CMD for the outer 
regions of the cluster. By eliminating these blends, we find that significant dusty mass loss occurs 
only at or near the tRGB. 



Two previous mid-IR studies of GCs detected several very red, faint sources dBoyer et all200^ 



Ita et al. 2007). Both studies offer the possibility that these sources are dusty, mass-losing stars. 



While the nature of the sources has yet to be established, we find that after weeding out pesky 
field stars and background galaxies via proper motion information, very few of the confirmed co Cen 
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cluster members reside in this region of the CMD (Fig. |9|, and significant dusty mass loss occurs 



only near the very tip of the RGB or AGB. In the case of NGC 362 (Ita et al. 2007), these faint 



red sources may be background AGB stars belonging to the Small Magellanic Cloud, however in 
other GCs, these sources are more likely to be background galaxies. 

If the current population in oj Cen is representative of the cluster's population at all times, then 
the cluster has lost at least 1-2 Mq of material (1 - 2 x 10"'* Mq of dust) through stellar mass 



loss since its last Galactic plane crossing, which occurred more than 3.4 x 10^ years ago dTayler &; 



Wood 19751. van Loon et al. (20071 find a few post- AGB stars in wCen. Assuming a post- AGB 



lifetime of several 10^ yr and 0.1 Mq of mass shed on the AGB, this could have produced up to 
100 Mq of lost mass. This comparison suggests that either these AGB stars lose on average less 
than 0.1 Mq on the AGB, or that our mass-loss rate estimates are too low (for instance because 
the wind velocity is actually higher and/or the dust-to-gas ratio is lower than assumed here). 



4.4. Background Galaxies 



Background galaxies (green points in Fig. [5]) dominate the faint, red plume in the [3.6] - [8] 



CMD due to strong PAH emission at 8 //m (Blum et al. 2006 Bolatto et al. 2007 1 . The coverage 
overlap between the 3.6 ^m and 8 //m mosaics is ~0.3 deg^, yielding a density of 2823 galaxies/deg^ . 
This number increases to 3174 galaxies/deg^ if we consider undetected galaxies behind the cluster 
center (an area of ~11 arcmin^, see Fig. [tJ. The red and faint region of the [3.6] - [24] CMD is 
also heavily populated with sources. If these 24 //m sources are all background galaxies, the galaxy 
density derived from this CMD is 2681 galaxies/deg'^ , which increases to 3015 galaxies/deg^ when 
considering sources located behind the cluster core. When including all sources with very red colors 
and faint magnitudes at all wavelengths, we find a total background galaxy density of ~1 x 10^ 
galaxies/deg^ (Fig. [tJ. 

To determine the expected background galaxy density, we queried the SWIRE extragalactic 



database (Lonsdale et al. 2004 1 for sources located in the same region of the CMD as our galaxy 
candidates. In the six SWIRE fields, there are «900 - 2200 sources/deg^ in the [3.6] - [8] CMD 
and 900 - 1900 sources/deg^ in the [3.6] - [24] CMD. Our images are slightly more sensitive than 
the SWIRE data, but when comparing identical regions of the CMD, the galaxy densities in our 
CMDs fall within the range of the SWIRE galaxy densities, strengthening the argument that the 
sources we have identified as galaxy candidates are indeed background galaxies. It is also possible 
that the faintest sources in the red [3.6] - [8] and [3.6] - [24] plumes are either stellar blends or 
regular RGB stars whose fluxes are poorly determined due to their proximity to the completion 
limit, causing them to have a non-zero color. 



4.5. 70 fim Sources 

Nine point-sources appear in the 70 fim mosaic, although none are the red, possibly dusty clus- 
ter members discussed in § 4.2 and § 4.3 One of these sources is resolved at the shorter wavelengths 
and appears to be a background spiral galaxy (2MASX J13272621-4746042; Fig.Il2]). This galaxy 
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shows a ring of PAH emission at 8 //m that is reminiscent of the famous star-forming "Ring of Fire" , 



seen in M31 (Barmby et al. 2006 Gordon et al. 2006). Two of 70 ^m sources are near confirmed 



cluster members ( van Leeuwen et al.||2000 van Loon et al.|2007D ; source 1 is 3.7" from LEID 27094, 
a distance that is within the half- width, half-maximum of the 70 /im PSF, and source 4 is 1.0" from 



LEID 46055. Optical spectra for both LEID stars are reported by van Loon et al. (2007) and reveal 
that both stars have [Fe/H] ~ -1.25 and T^s ^ 5000 K. LEID 27094 is CH strong and appears to 
be a normal RGB star, while LEID 46055 is CH and CN strong, and is Ba-rich, suggesting that it 
may be enriched in carbon from third dredge-up. Despite the proximity in position between these 
two LEID sources and sources 1 and 4, the LEID sources may be associated with warmer stars 
that appear in the IRAC mosaics as secondary stars blended with sources 1 and 4, and not with 
the corresponding 70 fim source. The remaining six 70 //m sources are not near any known cluster 
members. 

We measured 70 /xm fluxes for the eight point-sources with the aperture feature of Imexamine 
in version 2.12.2 of IRAF]^ Aperture radii were initially set to 3x the FWHM of each source and 
three iterations were performed to find the optimal radius. The median of the background was 
subtracted from the flux and the appropriate MIPS aperture correction was applied. Fluxes are 
listed in Table [7] and the 3.6 to 70 /im SEDs are shown in Figure 13 Aside from sources 1 and 4, 
all of the sources appear to be unaffected by blends at all IRAC wavelengths. 

For each 70 /xm source, we performed a reduced-x^ fit of a two-component blackbody curve 
that include 2MASS JHK fluxes to better define the warmer blackbody. Source 8 falls off of our 
IRAC images, so its IRAC fluxes were determined by using an earlier epoch of data from the Gehrz 
Guaranteed Time Observing Program with Spitzer (PID 132) with data reduction identical to that 
described in §[2j The fits take the photometric errors into account. Sources 3, 5, 7, and 8 are 
underconstrained for this fitting technique when we exclude the anomalously bright 5.8 and 8 /im 
points from the fit; for these sources, we decreased the number of fitting parameters by fixing the 
temperatures to 1000 to 14,000 K and 40 to 90 K for the blue and red blackbodies, respectively, and 
choosing the temperatures that gave the best reduced-x^- The IRAC and MIPS color corrections 
were applied by iterating the fit until the corrections remained constant within 1%. While a dual- 
blackbody fits well with most of the 70 fim sources, we point out that the dust around these sources 
may be emitting over a range of temperatures. 

Better fits were obtained by excluding the 8 /xm point, which is above the blue blackbody for 
all eight sources. To a lesser degree, the 5.8 fim point is also systematically brighter than the blue 
blackbody component, suggesting that there is a range of dust temperatures causing warm dust 
to radiate above the blackbody at 5.8 /xm and perhaps also at 8 fim. However, the 5.8 fiui points 
are still included in the fit for source 3 in order to keep the fit properly constrained. The resulting 



blackbodies are overplotted in Figure 13 and fitting results are shown in Table ^ 



The emission features at 5.8 and 8 fim, where PAH emission dominates, and the relatively flat 
SEDs of the non-blended sources suggests that all or most of these sources are background galaxies 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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(for examples of Spitzer galaxy SEDs, see Dale et al.|2005 ). For comparison, an SED of the resolved 



spiral galaxy is shown in Figure 14 We determined the flux by placing square apertures of roughly 
45" X 50" on the galaxy and subtracting the median background. As with source 8, the galaxy falls 
off of our IRAC coverage, so data from PID 132 is used to determine its IRAC fluxes. This SED 
is similar to the other eight, except that it lacks the slight excess emission at 5.8 /im that many of 
the other sources exhibit. 

The cold dust component of the blackbody fits for the eight unidentified sources range in tem- 
perature from ~50 K to 68 K - a temperature commensurate with these objects being post-thermal 
pulse stars, stars suffering from episodic mass loss on the RGB, or post-AGB stars. However, the 
temperatures of all eight sources are also consistent with dust in star forming regions that dominate 
the galactic far-IR emission. At a temperature of 48.8 it 0.8 K, the cold blackbody component 
of the resolved galaxy is colder than the other eight 70 iim sources, which may suggest that the 
warmer sources are of a different nature. Although, with so few constraints in the blackbody fits, 
one might argue that the temperatures of all nine sources are consistent with each other. 

If the 70 /xm sources are dusty AGE stars, the flux at 8 /um may be included in the wing of 
a broad emission feature centered at 9.5 ^m due to amorphous silicate grains that overlaps with 
the IRAC 8 ^m bandpass. The 9.5 /im feature is also detected in the 8-13 nm spectrum of 



a particularly red, large- amplitude variable star in 47 Tuc (van Loon et al. 2006a; Lebzelter et 



al. 2006). The positions of the 70 /um sources on the [8] - [24] CMD are marked in Figure [8p. 
All sources lie redward of the RGB and near the brighter end of the CMD. However, the 3.6 
and 4.5 //m magnitudes of sources 2, 3, 6, 7, and 8 (triangles) are fainter than 12.5 magnitudes, 
suggesting that these sources are indeed galaxies. Sources 1 and 4 lie on the brighter end of the 
RGB, but they are blended with the aforementioned LEID sources at the shorter wavelengths, 
which will cause an artificial increase in brightness. While near-IR redshift measurements or radial 
velocity measurements would conclusively determine the nature of these sources, analysis of the 
Spitzer data strongly suggests that they are all background galaxies. 



4.6. The Intracluster Medium 



Despite ongoing stellar mass loss in oj Cen, there are no known detections of any phase of ICM 



material in the cluster. The wind velocities discussed in §4.3 are <C 10 km s~^ for metal-poor, 
dust-driven winds from low-mass AGB stars. Such velocities are much less than the escape velocity 
of Lij Cen at 1 rhm ("Wesc = 44 km s~"^; Gnedin et al. 2002), but if the dust and gas are not well 



coupled, then the dust velocity may be higher. RGB and AGB wind velocities measured from Ha 
emission wings can be as high as the cluster escape velocity, however it is likely that Ha wings are 
due to pulsation and/or chromospheric activity instead of the circumstellar shells of stars. Core 
velocities originate in the chromospheres and are typically lower (< 15 km s~^), but there have 



been at least two observations of high velocities from the IR He I line (Dupree et al. 1992 Smith 



et al. 2004). 



If the winds are indeed slow and the total mass-loss rate discussed in S 4.3 has been constant 



since the last Galactic plane crossing, then, barring any ICM removal mechanisms and assuming 
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that all stars belong to the population with the lowest metallicity, 1-2 Mq of material should 
have accumulated in the cluster. This mass estimate includes only dusty mass loss on the RGB 
and AGB; mass loss from non-dusty stars on the RGB and from stars that lose significant mass in 
such a brief period of time that it may easily have been missed (e.g., associated with the formation 
of a planetary nebula) may cause this estimate to increase by a factor of two or three. Based on 
the total amount of AGB mass loss, which is the only source of dust, we expect a dust mass of Afd 
= Mgas ~ 1 ~ 2 X 10~^ Mq of dust. Since the estimate of the time since the last Galactic plane 
crossing is a lower hmit (> 3.4 x 10^ yr; [Tayler Wood||l975[ ), the predicted ICM dust mass may 
be higher. 



The average 1 a sensitivity of our 24 /im mosaic is ~2.5 x 10 MJy/str. Following Evans et 



al. (20031, the relationship between the dust mass and the expected IR flux density is: 



4.79 X 10"^^ /i.(mJy) 



kpc 



(3) 



where I?kpc is the distance to wCen, k^, is the dust absorption coefficient (cm g^ ), is the 
dust temperature (K), and B{v,Tci) is the Planck function in cgs units. Assuming a standard 
Mathis-Rumpl-Nordsieck dust distribution ( Mathis et al. [11977 1 and an ISM- type dust composition 
of graphite and silicate grains, is taken to be 56 it 11 cm^ g~^ (Ossenkopf &: Henning 1994 1, 
although this value could be up to a decade larger. A cloud with 1 - 2 x 10~^ Mq of dust and 
an equilibrium temperature of 70 K would yield a 24 ^m flux density of 23 - 45 mJy. With our 
sensitivity, we are capable of detecting 23 mJy to the 3 a level in our 24 ^m mosaic if all of the 
dust is uniformly distributed within an area of 2' x 2', or a radius of ~ 1/4 rhm A small, centrally 
located cloud at the 3 a level would be detectable over the background stellar emission in the 24 



and 70 //m mosaics. However, unlike in M15 (Boyer et al. 20061, our MIPS images do not show a 
large dust cloud located at the center of uj Gen. 

While there is no central IGM dust cloud in a; Gen, there is faint, diffuse emission covering 
nearly the entire 24 ^m image to the south of the cluster center, including several small regions 
of concentrated extended emission on the outskirts of the cluster that appear near the 3 a limit 



(Figs. 15 and 16). Feature A (Fig. 17) appears at 8, 24, and 70 /um, features B and G appear at 8 



and 24 ^um, and feature D appears only at 24 /im, as the coverage is not identical across all bands. 
The probability that these features are artefacts is quite low as they each appear in more than 95% 
of the BGD images at various array positions. 

It is uncertain whether features A - D are associated with the cluster or are part of the 
foreground. The morphology of feature A provides the strongest evidence for cluster association, 
as the warmer 8 /im emission lies closer to the cluster center than the cooler 24 and 70 ^m. emission 



(Figs. 16 and 17). Moreover, the warmer portion is approximately perpendicular to the direction 



towards the cluster center. This configuration suggests the possibilities that either a warm wind is 
pushing the IGM out of the cluster or that the 8 /xm arc is optically thick and is being irradiated 
by the cluster's stellar radiation field with cold dust in its shadow. A determination of the dust 
temperature would aid in ascertaining the origin of feature A. Unfortunately, only very small 



portions of feature A are present in the Spitzer mosaics at all wavebands (Fig. 17), and the limited 
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coverage makes it impossible to reliably determine the temperature of the dust from an SED. 



The case for an ICM origin for features A - D is strengthened by our discovery of an unreported 



H I cloud near the velocity of uj Cen (Fig. 18) in public data from the H I Parkes All-Sky Survey 



(HIPASS; [Barnes et alT][200T| ), which may have also been detected by [Smith et al.l ([1990[) in earlier 
observations. The center of this H I feature is not spatially coincident with any of the dust features, 
but rather lies «15' to the southeast of the cluster center, just on the edge of the coverage of our 
MIPS mosaic (Fig. 15). The velocity of the H I feature (i^LSR ~ 190 km s~^) suggests a connection 
230 km s^^l 



Smith et al. (1990) argue instead that, because the feature is 



to a; Cen (ulsr 

detected on only one side of the cluster, it is more likely to be a small high-velocity cloud that 

200 km s^^). While an association 



is part of the Magellanic Stream (Bajaja et al. 



1985 



between the H I feature and the Magellanic Stream is a possibility, similar offset H I clouds have 
been discovered in larger dSph galaxies (Bouchard et al. 2003 2006 Young et al.[[2006 ), so finding 
a cloud that is associated with oj Cen may not come as a surprise. In M15, dust and H I trace each 



other in the ICM (van Loon et al. 2006b). Therefore, if the H I feature is indeed associated with 



Lo Cen, it is reasonable to assume that dust is also present. Since both the dust and the H I feature 
are located well outside 1 rhnu it is possible that we are seeing the ICM in the process of escaping 
the cluster. 



5. Conclusions 

Our Spitzer multi-wavelength and dual-epoch study of the most massive Galactic globular 
cluster, a; Cen, and the resulting point-source catalog provide the most complete mid- to far-IR 
atlas of any GC to date. Despite lo Gen's rather unusual stellar population, the Spitzer CMDs show 
little structure beyond the division between the RGB stars belonging to the cluster and background 
galaxies. The HB stars and RGB-a stars are not separated from the RGB, except for a slightly 
blue [3.6] - [4.5] color for the RGB-a. 

Nine sources are detected at 70 fim, one of which is a resolved spiral galaxy, and two of which 
are spatially coincident with optical sources that are known cluster members. If emission at 5.8 and 
8 /um in the SEDs of the eight point sources is PAH emission, then these sources are background 
galaxies. 



By cross-referencing with a catalog of optical spectra (van Loon et al. 2007), we confirm that 
super-Ba-rich stars delineate the tip of the AGB, likely due to a third dredge-up, and find that three 
M-type stars not only have strong red excess, but are brighter than the 24 /im tRGB. Attributing 
the excess IR luminosity to dust suggests that the three bright M-type stars dominate the cluster's 
cumulative mass-loss rate (gas -|- dust) of 2.9 - 4.2 x 10"'^ Mq yi"^. If the cluster mass-loss rate 
has remained constant since the last Galactic plane crossing, then the cluster has lost at least 1 - 
2 Mq of dust-traced material over the last 3.4 x 10® years. 

The dusty stars in our sample range over ~1.5 dex in metallicity, suggesting that dust produc- 
tion is not inhibited even in stars with very low metal contents. In addition, our results show that, 
in LO Cen, significant dusty mass loss occurs near the very tip of the RGB or AGB, concentrated in 
only a few individual stars, with little evidence for such mass loss lower on the RGB. 
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While our images are sensitive enough to detect the amount of dust expected in uj Cen based on 

the current mass-loss rate, no obvious ICM dust clouds are apparent in the cluster center at 24 or 
70 /xm. We do find several regions of faint extended 24 /im emission on the outskirts of the cluster, 
but it is unclear whether these clouds are associated with the ICM or are in the foreground. If the 
dust features are associated with the cluster, they are well outside the cluster half-mass radius, and 
could therefore be ICM material that is in the process of escaping the cluster. 
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Fig. 1. — Three-color images of uj Cen. Upper panel: Blue is 3.6 fim, green is 4.5 fim and red is 
5.8 fim. Lower Panel: Blue is 3.6 fxm, green is 8 /xm and red is 24 fim. The black "X" marks the 
cluster center. Dusty stars and/or background galaxies become visible at 8 and 24 fim and appear 
red in the lower panel. No obvious intracluster medium dust clouds are visible. 
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Fig. 2. — Photometric errors. The x-axis is binned by 0.1 magnitudes, and the y-axis is binned 
by 0.2 percentage points. The dashed hnes mark the Hmit where photometry is 50% complete (see 
Table. [2]). The errors for the bulk of the stellar population are low, but climb steadily at the faintest 
magnitudes. 
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Fig. 3. — Color-color diagram showing which individual sources are affected by blending at 24 fim. 
We convolved the 8 /im image with the 24 ^m PSF and extracted point sources from the resulting 
image. Each panel shows the entire population in small points, with the 8 /um magnitude range 
indicated in each panel plotted in large, filled circles. Sources that have red [8] - [convolved 8] 
colors are affected by blending, whereas sources with red [8] - [24] colors and [8] - [convolved 8] ^ 
are truly red and likely harbor dust. We find that up to 25% of the population with [8] - [24] > 
0.1 is affected by blending. 




Fig. 4. — Luminosity functions for 3.6, 4.5, 5.8, 8, and 24 /im. The bin size is 0.2 magnitudes 
and fluxes are corrected for reddening {E{B — V) = 0.11; Lub 2002 1. In the upper panels, only 
sources detected in at least two bands are plotted and known non-members are excluded. In the 
lower panels, only known cluster members are plotted ( van Leeuwen et aL]|2000 ) . The dashed lines 
show the 3 a sensitivity limit for each band (at 3.6 and 4.5 /^m, the sensitivity limit is fainter than 
the plot area shown), and the dotted lines mark the 50% completeness limits. The photometry 
is nearly complete down to the 3 a sensitivity limit at 5.8 - 24 /xm, but is less complete at the 
shortest wavelengths due to severe stellar crowding in the cluster core. In the upper panels, the 
arrows indicate changes of slope that show the position of the RGB Bump/Clump and the HB. In 
the lower panels, the tRGB is marked as the point where the source counts drop by a factor of ~4. 




Fig. 5. — CMDs for oj Cen. Fluxes are corrected for reddening {E{B — F) = 0.11; Lub 2002). B 
magnitudes are from van Leeuwen et al. (20001 and J magnitudes are from Sollima et al. (2004 1 . 
The solid lines mark where the photometry is 50% complete in the case of Spitzer or mark the 
detection limit from van Leeuwen et al. (20001 and Sollima et al. (2004). The brightest confirmed 
cluster members at IRAC wavelengths are marked with large filled circles in each panel to show 
their evolution over color-magnitude space (non-members remain small dots). The three most 
extreme M-type stars are labeled with their LEID numbers. 
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Fig. 6. — CMDs for uj Cen. Same as Figure |5] with different source types plotted in colors. Black 
points are confirmed cluster members (van Leeuwen et al. 2000), gray points are not included in 
the van Leeuwen et al. (20001 proper motion (PM) catalog, green points are possible background 
galaxies, blue points are HB stars, and magenta points are sources belonging to the RGB-a. HB 
and RGB-a stars were identified by their positions on the B - [3.6] CMD. At Spitzer wavelengths, 
the HB does not separate from the RGB as in optical CMDs, nor does the RGB-a except for a 
trend towards slightly blue [3.6] - [4.5] colors. Background galaxies are identified by their [3.6] ~ 
[8] and [8] - [24] colors and appear very red at 8 /xm due to PAH emission and at 24 /xm due to 
dust emission. 
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Radius (arcmin) 

Fig. 7. — Source density profiles for uo Cen. Tlie tliick line shows possible HB stars in the top 
panel, AGB stars in the middle panel, and galaxy candidates in the lower panel. For comparison, 
the density profile for all cluster members is shifted down to match each profile at r = 8.5' and 
plotted in arbitrary units (a.u.) as a thin line. The AGB stars appear to be slightly more centrally 
concentrated than the general uj Cen population, and there is a paucity of HB stars and galaxy 
candidates in the cluster center. The HB and AGB profiles might be explained by the lower and 
higher flux densities of HB and AGB stars, respectively, than the general stellar population rather 
than a real difference between the profiles. 
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Fig. 8. — [8] - [24] CMDs showing the location of various stellar populations. Panel (a) shows 
super-Ba-rich and CN-rich stars as identified by their optical spectra ( van Loon et al.||2007 |; super- 
Ba-rich stars delineate the upper part of the AGB possibly due to third dredge-up. Panel (b) shows 
the 70 sources. Two of these are spatially coincident with cluster members (LEID sources), 
and one is a resolved spiral galaxy (2MASX J13272621-4 746042). The numbers to the right of 
each symbol correspond to the source numbers in Tables [7] and [8j All sources show excess 24 fj,m 
emission. 
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Fig. 9. — [8] - [24] CMDs for uj Cen. Cluster members are plotted in the left panel, non-members 
are plotted in the lower right panel, and sources lacking proper motion measurements are plotted 



in the upper right panel. Sources that are potentially blended at 24 fim (see §3.2) are not plotted. 



Regions including AGB mass-losing stars and dusty RGB stars are separated by dashed lines in 
the left panel. The colors corresponding to two sample optical depths (ri and T2) are also marked 
on the left. See Table [5] for the values of r for different dust compositions. M-type stars, carbon 
stars, and post-AGB stars identified by their optical spectra ( van Loon et al.|2007 | are also labeled, 
and the most extreme examples are also labeled with their LEID numbers. Almost all bright and 
red cluster members in this CMD have been identified against field stars, with only one bright star 
remaining for which no proper motion information is available. This source may be LEID 55017, 
which has a large proper motion and is not a cluster member. The bright luminosities and red colors 
of the three brightest M-type stars suggest that these three stars dominate the cluster mass-loss 
rate. All stars to the right of the vertical dashed line are included in the cumulative mass-loss rate 

(Fig.fnf. 
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Fig. 10. — Metallicity versus average [8] - [24] color. The color is a weighted average of all stars 
with a particular metallicity. Metallicities are determined by optical spectroscopy (van Loon et al. 
2007). Dust production does not appear to be inhibited at low metallicity. 
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Fig. 11. — Cumulative mass- loss rate (S M) for uj Cen. The upper panel shows the ratio of dusty to 
non-dusty stars (/d)- The dotted line includes all observed members, while the solid line excludes 
potential source blends. The bin size is 0.5 magnitudes, but the large bin between 4.5 < [24] < 7 
includes only the 3 brightest M-type stars, all of which are dusty. In the lower panel, individual 



mass-loss rates are derived from Groenewegen (20061 models scaled to reflect the luminosity and 



metallicity dependence of the dust-to-gas ratio, and consequently, the wind speed. Sources to the 
right of the vertical dashed line in Figure |9] are included in the summation and plotted with a solid 
line. The mass-loss rates are computed for four dust composition/temperature combinations, the 
most extreme of which are plotted; colder silicate and warmer AlOx grains fall intermediate to the 
plotted lines. The dashed line shows the cumulative mass-loss rate obtained if potential blends are 
included, which can increase the total mass-loss rate by 12% - 17% even in the relatively sparse 
environment of co Cen compared to other GCs. The sharp increases in S M beyond [24] ~ 10.5 are 
due to a few extremely red sources ([8] - [24] > 2). The total mass-loss rate for the cluster is 2.9 
- 4.2 X 10~'^ Mq yr~-^, 63% - 66% of which comes solely from the three brightest M-type stars, 
depending on the dust temperature and composition. 
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Fig. 12.— Three-color images of a resolved spiral galaxy (2MASX J13272621-4 746042). Left: Blue 
is 3.6 /im, green is 4.5 //m, and red is 8 nm. Right: Blue is 3.6 //m, green is 8 /im, and red is 24 /im. 
8 fim PAH emission appears red in the left panel and green in the right panel. Dust appears red in 
the right panel. The SED of this galaxy (Fig. 14 ) resembles the SEDs of the other 70 fim sources 
(Fig.fll. 
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Fig. 13. — SEDs of 70 sources. Fluxes are corrected for reddening {E{B — V) = 0.11; Lub 



2002). Sources 1 and 4 are spatially coincident with known cluster members and are potentially 
affected by blending between the known source and a fainter source that is the real counterpart of 
the 70 /xm emission. The 5.8 and/or 8 fiui emission in sources 2, 3, 6, 7, and 8 lies significantly 
above the blackbody and is excluded from the fit. Fitting results are summarized in Table |8j While 
the cold blackbody temperatures are consistent with a detached dust envelope such as those seen 
around post-AGB stars, they are also typical of the dust in star forming regions that dominate the 
galactic far-IR emission. The strong excess at 8 fini and the slight excess at 5.8 /xm could be due 
to either silicate emission or PAH emission. 
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Fig. 14.— SED of galaxy 2MASX J13272621-4746042. The 8 excess is suggestive of PAH 
emission. Blackbody fitting results are summarized in Table [8| This SED resembles the SEDs of 
the other eight 70 /im sources, suggesting that they are also galaxies. 
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Fig. 15. — 24 fim image showing the locations of four regions of extended emission. The boxed 
regions are shown in Figure 16 overlain with contours. Region A is shown in more detail in Figure 17 



An H 1 feature is marked with an "X" (Fig. 18). Low-level diffuse emission is also visible covering 
most of the image south of the cluster center. All dust clouds are outside of the cluster half-mass 
radius, suggesting that if they are part of the ICM, they may be in the process of leaving the 
cluster. 
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Fig. 16. — Residual 24 fim image showing several regions of extended emission. Contours smoothed 
with a 25 pixel {^60") window are overlaid to show the morphology of each region. The bright 
region on the edge of the frame in panel C is a data artefact. 
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Fig. 17. — (a) Residual 8 /im mosaic, (b) Residual 24 fim mosaic, (c) 70 /xm mosaic. Each panel 
shows a close-up of possible ICM feature A located ~20' from the cluster center. The arrows point 
towards the center of the cluster. The warm 8 ^m feature lies closer to the cluster center than 
the cooler feature at 70 iim. Both features appear at 24 /im. This morphology suggests a possible 
association with the cluster. Note that the region of bright emission to the lower left of feature A 
at 8 jum is strongly affected by residual emission from an extremely bright source, and is not likely 
a real extended feature. 
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Fig. 18. — H I detection near the velocity of uj Cen (fLSR ~ 230 km s ) in the HIPASS survey. 



This feature is located ~15' southeast of the cluster center (see Fig. 15). The presence of H I 
suggests that dust might also be present in the ICM of the cluster. 
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Table 1. Observational Summary 



Parameter 


IRAC 


MIPS 


A (^im) 


3.6, 4.5, 5.8 and 8 


24 and 70 


Program ID 


20648 


20648 


Date (UT) 


2006 Mar 26 


2006 Feb 22 and Mar 02 


AORkey 


14955264 


(1495)5520, 5776, 6032, 6288, 






6544, and 6800 


Coverage 


35' X 35' 


24 ;um: 35' X 55' 






70 fj,m: 21' X 55' 


Single Frame Exposure Time (s) 


10.4 


24 ^m: 9.96 






70 Atm: 10.49 


Average Depth (no. frames) 


10 


24 fim: 67 






70 urn: 36 


Pixel Size (arcsec pixel" ^) 


1.22 


24 nm: 1.25 X 1.30 






70 tJ.m: 4.93 x 5.03 


PSF size 


3.6 fj,m: 1.7" x 1.6" 


24 nm: 5.3" X 5.9" 




4.5 Atm: 1.6" X 1.6" 






5.8 fim: 1.6" X 1.7" 






<s Jim: l.<s" X 1,8" 





Note. — Observations are centered at R.A. = 13^26™58!4, decl. = -47°45'21'.'8 [J2000.0] 
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Table 2. Completeness Limits for Different Cluster Radii 



Wavelength 


50% Complete 


90% Complete 




(mag) 


(mag) 


3.6 


16.7 ± 0.1 


15.9 ± 0.1 


4.5 


16.7 ± 0.1 


15.9 ± 0.1 


5.8 


16.0 ± 0.1 


14.5 ± 0.1 


8 


15.8 ± 0.1 


14.3 ± 0.1 


24 


11.9 ± 0.1 


11.6 ± 0.1 



Note. — Completeness limits were determined 
from false star tests. These limits reflect the aver- 
age over the entire cluster coverage; photometry 
is less complete in the core than on the outskirts 
of the cluster. 
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Table 3. Sample Source Catalog 



No. 


Source ID'^ 




App 


arcnt Magn 


tudc 




Member- 


Cand- 


Note'' 


LEID" 






3.6 /xm 


4.5 


5.8 /im 


8 


24 ^^m 


ship (%)'' 


idate'^ 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


18562 


J132626. 28-472520. 06 


8.49(05) 


8.52(06) 


8.53(04) 


8.47(04) 


8.48(05) 


100 




CN 


35090 


18563 


J132626. 28-473229. 59 


14.04(09) 


14.16(07) 


14.24(07) 


14.32(07) 












18564 


J132626. 28-471144. 58 


16.58(08) 


16.37(09) 


15.86(13) 














18565 


J132626. 28-474102. 18 


16.77(09) 


16.68(16) 


16.33(18) 














18566 


J132626. 29-473737.68 


17.51(12) 


17.71(46) 
















18567 


J132626. 29-471200. 23 


17.20(11) 


17.13(20) 
















18568 


J132626. 29-473049. 11 


13.59(05) 


13.58(07) 


13.68(06) 






100 


HB 




47097 


18569 


J132626. 29-472651. 06 


13.23(05) 


13.20(06) 


13.18(05) 


13.16(10) 




98 






39111 


18570 


J132626. 30-471915. 97 


17.04(09) 


17.00(22) 
















18571 


J132626. 30-473515. 89 


17.30(14) 


17.55(27) 
















18572 


J132626. 30-471416. 24 


16.00(06) 




16.21(19) 














18573 


J132626. 30-474159. 18 


13.22(06) 


13.24(06) 


13.28(04) 






86 




CN 


72019 


18574 


J132626. 31 -471627.20 


9.34(04) 


9.28(06) 


9.31(05) 


9.29(04) 


9.30(06) 


IS 




P 


16018 


18575 


J132626. 31-473111. 71 


14.94(08) 


15.06(11) 
















18576 


J132626. 31 -472320.08 


13.98(06) 


13.99(05) 


14.00(05) 


14.02(06) 












18577 


J132626.31 -474246. 15 


17.80(20) 


17.95(52) 




16.08(35) 






Galaxy 






18578 


J132626. 32-473543. 47 


13.12(04) 


13.18(09) 


13.11(08) 


13.18(08) 












18579 


J132626. 32-472534.32 


10.13(05) 


10.15(06) 


10.15(05) 


10.16(04) 


10.18(09) 


100 






36061 


18580 


J132626. 32-474718. 12 




15.72(08) 




15.67(26) 












18581 


J132626. 33-473132. 57 


12.23(08) 


12.28(07) 


12.24(06) 


12.02(06) 




100 


RGB-a 




49082 



Note. — Magnitude uncertainties arc quoted in parentheses. For example, a magnitude of 8.49 zb 0.05 is listed as 8.49(05). 

^ To save space, the Source ID is truncated to show only the source coordinates. The full electronic catalog, available through CDS, includes 
the Spitzer pretix (SSTOCEN). 

^ Membership likelihoods arc quoted as a percentage from jvan Leeuwen et al.| l |2000[ . 

^ Candidate AGB, RGB-a, and HB stars and candidate galaxies arc identified by their positions on optical and Spitzer CMDs. 

The Note column specifics whether a source is CN- or super-Ba-rich and/or is a M-type (M), carbon (C), or post- AGB (P) star, as indicated 
by optical spectra Jvan Loon et al. |2007t . The note also indicates if a source is potentially affected by blending at 24 fim (B24, see § |3.2[ l . 

° The LEID identification numbers are assigned to Spitzer sources via cross-identification with[van Leeuwen et al.H2000|. 



Table 4. Spitzer MAGNITUDES OF M-Type, Carbon, and post-AGB Stars identified by 

VAN Loon et al.| (|2007[) 



LEID 


Source ID^ 


m3.6 


"14.5 


ms.s 


mg 


m24 


Source Type 



14043 


J132750 


37- 


471542 


79 


11 


75(06) 






11 


77(04) 










Carbon Star 


16018*^ 


J132626 


31- 


471627 


20 


9 


34(04) 


9 


28(06) 


9 


31(05) 


9 


29(04) 


9 


30(06) 


Post-AGB 


30020 


J132549 


95- 


472259 


76 


12 


62(05) 


12 


62(06) 


12 


61(06) 


12 


66(04) 






Post-AGB 


32015 


J132538 


44- 


472401 


79 


10 


26(05) 


10 


23(07) 


10 


22(07) 


10 


27(04) 


10 


27(09) 


Post-AGB 


32029 


J132605 


17- 


472342 


43 


9 


11(05) 


8 


96(05) 


8 


82(04) 


8 


81(04) 


8 


48(05) 


Post-AGB 


32059 


J132628 


13- 


472340 


56 


8 


45(05) 


8 


52(06) 


8 


52(04) 


8 


40(05) 


8 


40(05) 


Carbon Star 


33062 


J132630 


19- 


472427 


89 


6 


96(05) 


6 


88(06) 


6 


74(05) 


6 


27(04) 


4 


98(05) 


M-Type 


35094 


J132628 


81- 


472523 


60 


8 


43(05) 


8 


56(07) 


8 


48(05) 


8 


50(05) 


8 


43(05) 


M-Type 


35250 


J132737 


72- 


472517 


37 


7 


39(05) 


7 


45(06) 


7 


35(05) 


7 


23(04) 


6 


61(05) 


M-Type 


35252 


J132738 


29- 


-472505 


60 


14 


21(04) 


14 


07(10) 


14 


49(08) 


14 


15(10) 






M-Type 


41071 


J132614 


42- 


472805 


39 


9 


82(04) 


9 


93(06) 


9 


89(05) 


9 


81(04) 


9 


97(08) 


Carbon Star 


42044 


J132605 


35- 


472820 


80 


7 


59(04) 


7 


74(06) 


7 


67(04) 


7 


55(05) 


7 


33(05) 


M-Type 


43105 


J132627 


22- 


472847 


47 


9 


56(04) 


9 


54(06) 


9 


46(04) 


9 


44(04) 






Post-AGB 


44262 


J132646 


37- 


472930 


35 


7 


05(05) 


6 


85(05) 


6 


68(04) 


6 


42(04) 


5 


85(05) 


M-Type 


44484 


J132726 


37- 


472916 


93 


8 


44(06) 


8 


63(05) 


8 


53(05) 


8 


51(05) 


8 


44(05) 


M-Type 


52030 


J132601 


59- 


473306 


08 


7 


40(05) 


7 


48(06) 


7 


44(05) 


7 


32(04) 


7 


33(05) 


Carbon Star 


53019 


J132544 


03- 


473324 


79 


13 


03(05) 


13 


07(06) 


13 


03(04) 


12 


77(10) 






Carbon Star 



Note. — Magnitude uncertainties are quoted in parentheses. For example, a magnitude of 11.75 it 0.06 is 
listed as 11.75(06). 

To save space, the Source ID is truncated to show only the source coordinates. The full electronic catalog 
(available through CDS) includes the Spitzer prefix (SSTOCEN). 

b 



Fehrenbach's Star jPehrenbach fc Duflot|l962 L 
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Table 5. Sample Optical Depths Corresponding to [8] - [24] = 0.25 and 1.7 





Silicat 


e Dust 


AlO, 


Dust 


Optical Depth 


[8] - [24] T = 3297 K 


T = 3850 K 


T = 3297 K 


T = 3850 K 


n 


0.25 0.013 


0.001 


0.030 


0.021 


T2 


1.70 0.120 


0.056 


0.220 


0.090 



Note. — Ti and T2 are examples of optical depth from |Groenewegen ( 2006[ ) that 
correspond to [8] - [24] colors of 0.25 and 1.7. These colors are marked in Figure |9] 
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Table 6. Mass-Loss Rates of the Three Brightest M-type Stars 



M {Mq yr-i) 
Silicate Dust AlOx Dust 



LEID 


Source ID 


T = 3297 K 


T = 


3850 K 


T = 


3297 K 


T = 


3850 K 


33062 


SSTOCEN J132630. 19-472427.89 


1.7 X 10"^ 


1.3 


X 10-^ 


2.1 


X 10-^ 


1.5 


X 10-^ 


44262 


SSTOCEN J132646. 37-472930.35 


3.4 X 10-** 


2.8 


X 10-8 


4.4 


X 10-8 


3.5 


X 10-8 


35250 


SSTOCEN J132737.72-472517.37 


2.3 X 10-8 


1.8 


X 10-8 


2.8 


X 10-8 


2.2 


X 10-8 




Total M for all dusty members 


3.4 X 10"^ 


2.9 


X 10-^ 


4.2 


X 10-^ 


3.4 


X 10-^ 
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Table 7. Fluxes of 70 fiM Sources 













Flux 


(mjy) 






Source # 


Source ID 




3.6 


4.5 fJ.m 


5.8 


8 


24 ^im 


70 


1 


SSTOCEN J132723.86- 


-472130.00 


15.94(0.32) 


10.23(0.20) 


6.66(0.13) 


5.59(0.11) 


2.42(0.05) 


75.79(17.31) 


2 


SSTOCEN J132521.01- 


-472623.43 


1.96(0.04) 


1.49(0.03) 


1.96(0.04) 


12.33(0.25) 


14.33(0.29) 


226.70(8.18) 


3 


SSTOCEN J132523.57- 


-472813.61 




0.17(0.00) 


0.16(0.00) 


0.79(0.02) 


2.82(0.06) 


246.73(18.94) 


4 


SSTOCEN J132615.23- 


-473024.56 


46.46(0.93) 


30.02(0.60) 


19.52(0.39) 


11.56(0.23) 


10.74(0.47) 


158.02(7.96) 


5 


SSTOCEN J132426.26- 


-473905.88 










0.74(0.05) 


82.01(9.44) 


6 


SSTOCEN J132714.06- 


-474112.83 


0.93(0.02) 


0.63(0.01) 


0.60(0.01) 


2.60(0.05) 


3.54(0.07) 


53.04(11.19) 


7 


SSTOCEN J132655.87- 


-474445.02 


0.12(0.01) 


0.11(0.01) 


0.12(0.01) 


0.69(0.01) 


3.81(0.08) 


131.46(23.81) 


8 


SSTOCEN J132712.89- 


-474739.59 


0.07(0.01) 


0.05(0.01) 


0.05(0.01) 


0.06(0.01) 


0.75(0.02) 


107.50(23.54) 


ga 


2MASX J13272621- 


4746042 


2.92(0.07) 


3.38(0.09) 


3.06(0.08) 


9.27(0.21) 


5.20(0.21) 


617.08(73.84) 



Note. — Fluxes for 3.6 — 24 /^m were determined by PRF fitting with DAOphot and 70 fim fluxes were determined by aperture photometry 
with IRAF. Flux uncertainties are quoted in parentheses. For example, a flux of 15.94 zb 0.32 mJy is listed as 15.94(0.32). 

^ Source 9 is a resolved spiral galaxy (see Fig. |l2| , while sources 1 through 8 are point sources. 



Table 8. Fitting Results for 70 fiM Sources 



Source # 


Source ID 




xVdof 


Ti (K) 


T2 


(K) 


1 


SSTOCEN J132723.86 


-472130.00 


35.8 


4664 ± 


19 


58.5 


± 


1.7 


2 


SSTOCEN J132521.01 


-472623.43 


3.3 


2688 ± 


25 


67.5 


± 


0.5 


3 


SSTOCEN J132523.57 


-472813.61 


199.7 


1100 ± 


200^ 


55.0 


± 


10^ 


4 


SSTOCEN J132615.23 


-473024.56 


40.1 


4329 ± 


16 


66.5 


± 


0.8 


5 


SSTOCEN J132426.26 


-473905.88 


75.0 






50.0 


± 


10^ 


6 


SSTOCEN J132714.06 


-474112.83 


15.4 


2901 ± 


28 


68.0 


± 


2.0 


7 


SSTOCEN J132655.87 


-474445.02 


52.2 


1800 ± 


200^^ 


60.0 


± 


10^ 


8 


SSTOCEN J132712.89 


-474739.59 


15.7 


4000 ± 


200^ 


50.0 


± 


10^ 




2MASX J13272621- 


4746042 


2.5 


1084 ± 


29 


48.8 


± 


0.8 



^ The marked uncertainties represent twice the temperature increment in the range 
of fixed temperatures used for the underconstrained fits. 



^ Source 9 is a resolved spiral galaxy (see Fig. 12 1, while sources 1 through 8 are 
point sources. 



